Geometric magnetism in classical transport theory 



On 
On 



oo 



43 

o 

CD 



CO 
C3 



I 

C 

O 
O 



> 

m 

m 
o 

On 



o 



X 
S3 



Jochen Rau 

Max- Planck- Institut fur Physik komplexer Systeme, Bayreuther Strafie 40 Haus 16, 01187 Dresden, Germany 
(revised version, May 26, 1997; to appear in Phys. Rev. E) 



The effective dynamics of a slow classical system coupled to 
a fast chaotic environment is described by means of a master 
equation. We show how this approach permits a very simple 
derivation of geometric magnetism. 
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I. INTRODUCTION 

Consider a slow classical system S coupled, through its 
position, to a fast classical system F. If the fast motion 
is chaotic then F effectively acts as an "environment" 
which induces friction |]||] and also exerts other, non- 
dissipative reaction forces on the slow system S. As 
for the non-dissipative reaction, in the simplest "adi- 
abatic averaging" approximation [Q -the classical ana- 
logue of the Born-Oppenheimer approximation ||]— the 
fast motion's energy at given values of the slow coordi- 
nates serves as an external potential for the slow system; 
its gradient yields the "Born-Oppenheimer force." In 
the next approximation beyond this, there is a velocity- 
dependent correction which has the form of a magnetic 
force, and for which Berry and Robbins coined the name 
"geometric magnetism" [|6|. 

For the situation considered here, namely a fast chaotic 
environment, there are so far two alternative derivations 
of geometric magnetism, both due to Berry and Rob- 
bins: (i) by taking the (non- vanishing) classical limit of 
the corresponding quantum result H , where the appear- 
ance of geometric magnetism can be linked to the geo- 
metric phase |sj (whence its name); or (ii) in a purely 
classical context, by expanding the equation of motion 
for S around the Born-Oppenheimer limit in powers of 
the fast/slow time scale ratio, and identifying the first- 
order correction l|. 

The multiple-time-scale analysis |J used in the classi- 
cal derivation of geometric magnetism appears somewhat 
reminiscent of the old Chapman-Enskog method to derive 
transport equations Jh}. This is not a mere coincidence: 
After all, it should be possible to describe the evolution 
of any subsystem (here: S) coupled to an environment 
(here: F) by means of a transport ("master") equation; 
and, to be consistent, such a master equation should fea- 
ture a term representing geometric magnetism. It is the 
purpose of the present paper to show how, indeed, geo- 
metric magnetism arises naturally in a classical master 
equation for S. The derivation of geometric magnetism 



within this transport theory framework will turn out to 
be surprisingly simple. 

The paper is organized as follows. First 1 will sketch 
very briefly how one obtains transport equations by 
means of the Nakajima-Zwanzig projection technique 
|ll]-[l4| (Se c. |ll| ); for more details the reader is referred to 
textbooks ]15[ and a recent review pTq ] . Then this projec- 
tion technique is applied to the situation at hand, namely 
a slow system S coupled to a fast chaotic environment 
F, both taken to be classical (Sec. 111). The resulting 
master equation for S is generally non-Markovian, yet 
the separation of time scales and chaoticity of the fast 
motion permit us to take the Markovian limit. Along 
with the Born-Oppenheimer force, geometric magnetism 
then appears in a straightforward manner in the non- 
dissipative part of the effective slow dynamics. Finally, 
I shall conclude with a brief summary and several addi- 
tional remarks (Sec. 



[V) 



II. TRANSPORT EQUATIONS 

A powerful tool for the derivation of transport equa- 
tions is the Nakajima-Zwanzig projection technique 
pl| |Tg[ . Its main strength lies in the fact that by map- 
ping the influence of irrelevant degrees of freedom onto 
-among other features- a non-local behavior in time, it 
opens the way to the systematic exploitation of separated 
time scales and hence serves as a good starting point for 
powerful approximations such as the Markovian and qua- 
sistationary limits; furthermore, it permits one to discern 
easily the dissipative and non-dissipative parts of the ef- 
fective dynamics. 

When studying the dynamics of a complex system 
away from equilibrium one typically monitors the evo- 
lution of the expectation values 



g a (t) := ( P (t)\G a ) 



(1) 



of only a very small set of selected ( "relevant" ) observ- 
ables {G a }, which change according to 



g a (t)=i(p(t)\CG a ) 



(2) 



Here the meanings of p(t), G a and (|), as well as of the 
scalar product (; ) p which we shall use later, depend on 
whether the system under consideration is quantum or 
classical; they are summarized in Table I. The "Liouvil- 
lian" C takes the commutator with the Hamilton opera- 
tor H, 



iC = {%/%) [H, ■ 



(3) 
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for a quantum system, or generates a Lie dragging in the 
direction of the Hamiltonian vector Xh, 



with summation over a implied (Einstein convention), 
partition function 



iC = £; 



(4) 



Z(0) := (l|cxp[-A a (0)G Q ]) 



(9) 



for a classical system 17 , respectively. For simplicity 
we assume that the Hamiltonian and hence the Liouvil- 
lian, as well as the relevant observables, are not explicitly 
time-dependent . 

The right-hand side of the equation of motion (||) will 
generally depend not just on the selected, but also on 
all other ("irrelevant") degrees of freedom. In order to 
eliminate these and hence obtain a closed "transport 
equation" for the {g a (t)}, one employs a suitable projec- 
tion operator which projects arbitrary observables onto 
the subspace spanned by 1 and the relevant observables 
{G a }. The projector may depend on the current expec- 
tation values of the relevant observables and thus vary in 
time, V{t) = V[g a {t)], and is assumed to have the three 
properties (i) V(t) 2 = V{t); (ii) V(t)A = A if and only if 
A 6 span{l, G a }; and (iii) 



Pit) 



dP(t) 



dt 



A =0 V p(t),A 



(5) 



Its complement is denoted by Q(t) := 1 — V(t). One 
further defines an operator T(t', t) (again in the space of 
observables) by 



d_ 

dt 



7 T(t',i) = -iQ(t')CQ(t')T{t',t) 



(6) 



with initial condition T(t,t) = 1, which may be pictured 
as describing the evolution of the system's irrelevant de- 
grees of freedom. The equation of motion for the selected 
expectation values {g a (t)} can then be cast into the -still 
exact- form 

g a (t)=i(p(t)\T(t)CG a ) 



dt' {p(t')\V(t')CQ{t')T{t\ t)Q{t)CG a ) 



+ i(p(0)\Q(0)T(0,t)Q(t)CG a ) , (7) 

for any time t > 0. This constitutes the desired 
closed system of (possibly nonlinear) coupled integro- 
differential equations for the selected expectation values 
{g a (t)} if and only if (p(0)|Q(0) and with it the last 
("residual force") term vanishes. 

In many practical applications the initial state p(0) is 
not known exactly but characterized solely by the ini- 
tial expectation values {g a (0)} of the relevant observ- 
ables. From this insufficient information one generally 
constructs that distribution which maximizes the entropy 
S[p] := —k(p\ lnp) and hence can be considered "least bi- 
ased" or "maximally non-committal" with regard to the 
unmonitored degrees of freedom: It is the generalized 
canonical state 



and the Lagrange parameters {A a (0)} adjusted such as 
to yield the correct {g a (0}}- 

In an analogous fashion one defines a "relevant part" 



p iel (t) :=Z(t)- 1 exp[-\ a (t)G a ] 



(10) 



of the exact state p{t) at all times t, where p re \{0) = p(0) 
but generally p rc \(t) ^ p(t) for t > 0. There exists 
a unique time-dependent projector 7-"r(£), namely the 
projector orthogonal with respect to the time-dependent 
scalar product (; ) Prol (t), which has all required properties 
(i)-(iii) and, moreover, yields 



(p(t)\V R (t) = (p iel (t)\ 



(11) 



at all times. This special choice, originally proposed by 
Robertson @,^8| , has the important advantage that for 
initial states of the form (Q) it ensures (p(0)|Qr (0) = 
and so renders the transport equation (0) closed. We 
shall use the Robertson projector throughout the remain- 
der of the paper (and, for brevity, immediately drop the 
subscript 'R'). 

One principal feature of the transport equation (0) is 
that it is non-Mar kovian: Future expectation values of 
the selected observables are predicted on the basis of both 
their present values and their past history. There are two 
distinct time scales: (i) the scale T re i -or several scales 
{ T rei}~ 011 which the selected expectation values {g a (t)} 
evolve; and (ii) the "memory time" T mem which char- 
acterizes the length of the time interval that contributes 
significantly to the memory integral. Only if this memory 
time is small compared to the typical time scale on which 
the selected observables evolve, r mcm <ti t ic \, can mem- 
ory effects be neglected and predictions for the selected 
observables be based solely on their present values. One 
may then assume that in the memory term g a (t') ~ g a (t) 
and hence replace 



V[9a{t')\ 

(p(t')W(t') 

T(t',t) 



■P[9a(t)} , 

(p(t)\V(t) , 
exp[iQ(t)CQ(t) ■ (t - t') 



(12) 



(Markovian limit). Furthermore, at times t 3> r mcm it no 
longer matters for the dynamics when exactly the evolu- 
tion started, and hence in Eq. (^) the integration over 
the system's history may just as well extend from — oo to 
t, rather than from to t (quasistationary limit) M|. In 
the Markovian and quasistationary limits the equation of 
motion simplifies to 

g a (t) =i(p(t)\£ xel (t)G a ) 

-n(p(t)\V(t)CQ(t)6(Q(t)CQ(t))Q(t)CG a ) (13) 



p(0) = Z(0)- 1 exp[-A Q (0)G a ] 



(8) 



modulo residual force, where 
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£ iel (t) = V{t)CV{t) 



with 



+ - / dTV(t)£[C(t;r)-C(t;-T)]V(t) (14) 



(15) 



C{t;r) := exp[iQ(t)CQ(t)T] C 



pertaining to S, as well as the total energy which is a 
constant of the motion, are taken to be relevant; while 
the internal degrees of freedom of the environment and 
system-environment correlations are deemed irrelevant. 
This gives rise to a time-independent representation of 
the Robertson projector, 



denotes a -possibly time-dependent- effective Liouvillian 
for the relevant observables. 

Provided the evolution operator T is unitary with re- 



VA = 



1 



d E « 

=: (A)e 



J dzS(H - E)A 



spect to the scalar product (; ) Prol(t) then the first term in f or any observable A. Here d E := d/dE and 
the Markovian transport equation ( |l3| ) is non-dissipative. 
In this case dissipation stems entirely from the second 
term, which yields a non-negative entropy growth rate 



n := / dzO{H - E) ■ 



(21) 



(22) 



S\pt*{t)\ =k\ a (t)g a {t) 

= kir(Q£\ b G b ;5(Q£Q)Q£\ a G a ) Piol >0 (16) 

(_ff-theorem). 



III. EFFECTIVE FORCE CAUSED BY A FAST 
CHAOTIC ENVIRONMENT 

We now apply the above general results to a slow 
system S coupled to a fast chaotic, but not necessarily 
macroscopic, environment F. Both S and F are treated 
classically, and their state described in a phase space with 
canonical coordinates Z = {Q, P} pertaining to S and 
z = {q, p} pertaining to F, respectively. The full Hamil- 
ton function for the combined system S x F is taken to 
be of the form 



H(Z,z) = H s (Z)+h(Q,z) 



(17) 



where H§ governs the free dynamics of the system S, and 
h describes both the coupling -through the slow position 
Q- of S to the environment and the internal dynamics 
of the latter. Associated with the Hamilton function is a 
Liouvillian (|]) which we decompose 



£ Xu £ X n z 



into a part dragging along the slow coordinates 
dQ 



X 



H.Z 



(18) 



(19) 



with slow velocity V 1 = dHs/dPi and di := d/dQ l , and 
a part dragging along the fast coordinates, 



X 



H. 



E 



dh d dh d 
dpk dq k dq k dpi 



(20) 



At t = and hence, due to energy conservation, at 
all times the combined system S x F is assumed to 
have a sharp total energy E. For the purposes of the 
Nakajima-Zwanzig projection technique all observables 



its derivative <9eO may be interpreted as the surface of 
the microcanonical energy shell. 

The slow system's effective dynamics must be de- 
scribed with a transport equation of the form (Q) , which 
in general is non-Markovian and includes a residual force. 
Only if we take the initial state of the environment to be 
microcanonical, i.e., 



P(0)=as(0) x —6(H-E) 



(23) 



where ps(0) denotes the (arbitrary) initial state of S, 
then the residual force term vanishes . This assump- 
tion of a microcanonical distribution and the resultant 
omission of the residual force term amount to averaging 
over an entire ensemble of fast chaotic systems. However, 
even if the slow system is coupled to a single fast chaotic 
system the transport equation without residual force will 
describe the main global feature of the slow dynamics; 
the residual force causes only fluctuations around the av- 
erage trajectory ]21f| . 

The separation of time scales and chaoticity of the fast 
motion permit us to take the Markovian and quasista- 
tionary limits. Moreover, we focus on the non-dissipative 
part of the slow dynamics. The latter is governed by the 
effective Liouvillian (|l4|), which immediately yields the 
(Heiscnberg-picture) equation of motion 

G = {£x H ,z G)e 

+ 2 J dr ^ £x »^^ £x h,z(t) - £x h .z(-t)]G)e (24) 

for an arbitrary slow observable G(Z) p^ |. Here we have 
used £x H Z V = T ' £x H , = to replace £x H by £x H z . 
and defined the "rotated" Hamiltonian vector Xh,z{t) 
as in (jl^) but with components dragged along the fast 
coordinates: 



dih -> (dih) T := exp[r£ Xf , J (dih) 



(25) 



Upon choosing G — P we obtain an effective force, with 
components 



3 



1 f°° 

Pi = -{d l H) E V j J dT (djftdiQr ~ {dih)- T ]) E 



_ pBO + ^geo 



(26) 



The first term constitutes the usual Born-Oppenheimer 
force; while the second (integral) term gives rise to geo- 
metric magnetism: For an arbitrary function A(Z, z) it 
is 



(djA) E = dj(A) E + (d 3 H) E d E (A) E 



d E fl 



d E n-{AdjH) E 



(27) 



where djH := djH — (djH) E . This identity, together 
h, 

((dih) T ) E = ((dih)^ T ) E (28) 



with djH — djh. 



and 



yields 



{(d l h)_ T d ] h) E = (0 3 h) T d l h) E , (29) 



»« r ' ( 3 °) 



example for this is Jarzynski's force |2j| . Here, too, trans- 
port theory with its non-Markovian evolution equation 
(0) will furnish a good starting point for the systematic 
study of non-Markovian corrections. 

Finally, transport theory adds a somewhat new per- 
spective to Berry and Robbins' observation that geomet- 
ric magnetism is the "antisymmetric cousin of friction" 
|^| . In our formulation the memory term in Eq. (Q) gives 
rise, upon taking the Markovian and quasistationary lim- 
its, to both dissipative and non-dissipative parts of the 
dynamics; they appear as the parts symmetrized and an- 
tisymmetrized, resp., with respect to the history integra- 
tion variable r = (t — t 1 ). When evaluated for our model 
Hamiltonian and the special choice G = P these two 
parts translate into effective forces proportional to the 
slow velocity, one with a symmetric matrix of coefficients 
(friction: not considered in this paper), the other with an 
antisymmetric matrix (geometric magnetism: Eq. (31)). 
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with an antisymmetric matrix ("magnetic field") 
1 



Bi 



2d E n 



d f 



d E n I d.T((dih) T djh — {djh) T dih) E 
o 



(31) 



in agreement with the result by Berry and Robbins 



IV. CONCLUSION 

We have succeeded in deriving geometric magnetism 
within the framework of classical transport theory: Start- 
ing from the general formula (|l4|), the derivation turned 
out to be surprisingly simple. This means that treat- 
ing the fast chaotic system as an "environment" for the 
slow system and describing the dynamics of the latter 
with a master equation is a consistent and useful physi- 
cal picture. More generally, it shows that methods from 
transport theory are not limited to the description of 
macroscopic systems but apply as well to low-dimensional 
chaos. 

As a by-product we have obtained an equation of mo- 
tion for arbitrary slow observables (Eq. (p4|)) which is 
coordinate-free and applies to any form of the micro- 
scopic Hamilton function: It could serve as the starting 
point for interesting generalizations of the model Hamil- 
tonian considered here. Higher-order corrections to the 
Born-Oppenheimer force and geometric magnetism will 
presumably have to take into account memory effects; an 
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TABLE I. Various symbols used in transport theory 



generic symbol 


quantum 


classical 


state p 


statistical op. 


P 


phase space dist. p(n) 


observable G a 


Hermitian op. 


Ga 


real function G a (j^) 


(A\B) 


tr(i f B) 




J d% A* (w)B{ir) 


(A;B) P 






f d-K p(tt)A*(tt)B(-k) 



